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ABSTRACT 



Black holes grow by accreting matter from their surroundings. However, angular 
momentum provides an efficient natural barrier to accretion and so only the lowest 
angular momentum material will be available to feed the black holes. The standard 
sub-grid model for black hole accretion in galaxy formation simulations - based on 
the Bondi-Hoyle method - does not account for the angular momentum of accreting 
material, and so it is not clear how representative the black hole accretion rate esti- 
mated in this way is likely to be. In this paper we introduce a new sub-grid model for 
black hole accretion that naturally accounts for the angular momentum of accreting 
material. Both the black hole and its accretion disc are modelled as a composite ac- 
cretion disc particle. Gas particles are captured by the accretion disc particle if and 
only if their orbits bring them within its accretion radius R acc , at which point their 
mass is added to the accretion disc and feeds the black hole on a viscous timescale 
t v i SC . The resulting black hole accretion rate Mbh powers the accretion luminosity 
L acc oc Mbh, which drives black hole feedback. Using a series of controlled numeri- 
cal experiments, we demonstrate that our new accretion disc particle method is more 
physically self-consistent than the Bondi-Hoyle method. We also discuss the physical 
implications of the accretion disc particle method for systems with a high degree of 
rotational support, and we argue that the Mbh — cr relation in these systems should be 
offset from the relation for classical bulges and ellipticals, as appears to be observed. 



Key words: accretion: accretion discs 
methods: numerical 



galaxies: active - galaxies: formation 



1 INTRODUCTION 

Understanding how super- massive black holes at the centres 
of galaxies grow over cosmic time is one of the most im- 
portant yet challenging problems facing modellers of galaxy 
formation. Observationally there is clear and compelling 
evidence that in galaxies that host super-massive black 
holes the black hole mass Mbh correlates strongly with 
the stellar mass M* and velocity dispersion a of the host 
buke (e.g.|Magorrian et al.lll998l:lFerrarese fc Merrittll2000l: 



iGebhardt et all l2000l: iTremaine et al.1 120021 : lHaring fe Rixl 
120041 ; iGiiltekin et al.ll2009r ). Theoretically these correlations 
are widely interpreted as hallmarks of black hole feedback, 
which itself is a natural consequence of accretion onto the 



the feedback and the greater the impact on Mbh - Therefore, 
how one estimates Mbh is crucial because it governs not 
only the rate at which the black hole grows but also the 
strength of the black hole feedback. This is a particularly 
important problem because how black hole feeding and 
feedback is modelled can have a profo und impact on the 
predictions of how galaxies form (e.g. iBower et al. 1 120Q6I : 
ICroton et al]|2006h . 

The standard approach to estimating Mbh in 
gal axy formation simula tions is bas e d on the work 
of iBondi fe Hovkl (|l944h and iBondil dl952h (he r eafter 



black hole (e.g. ISilk fe Reed 1 19981 : lFabianlll999l : lKindl2Q03l : 



ISazonov et alTl2QQ5l : iKingl l20Q5h . In this picture, feedback 
acts to regulate the black hole's mass accretion rate Mbh 
by modifying the physical and dynamical state of gas in and 
around its host galaxy - so the greater Mbh, the stronger 
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the Bondi-Hoyle method; cf. iDi Matteo et al.1 120051 : 
ISoringel et aljEooBh . In the accretion problem as it was orig- 
inally formulated, a spherically symmetric accretion flow is 
captured gravitationally by a point-like accretor from a uni- 
form distribution of gas with zero angular momentum. Un- 
der these conditions, the accretion rate onto the accretor 
MBondi is proportional to the square of the black hole mass 
Mbh an d the gas density /?, and inversely proportional to 
the cube of the sound speed c s . This gives an accretion rate 
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MBondi oc M^p/c 3 s . The assumption in galaxy formation 
simulations i s that Mbh oc MB on d i (see, for example, the 
discussion in iBooth Schavell2009h . 

However, there are good physical reasons to believe 
that M Bon di cannot be representative of the true black 
hole accretio n rate M bh in an astrophysically realistic 
situation (cf. lKindl201Qh . First, the black hole is embedded 
in the gravitational potential of a galaxy that is orders 
of magnitude more massive than it; this means that the 
gravitational force acting on the accretion flow is dominated 
by the mass of the galaxy rather than the black hole and so 
MBondi will be a similar number of orders of magnitude off 
the true Mbh (we show this explicitly in Hobbs et al., in 
preparation). Second, any astrophysically realistic accretion 
now will have some angular momentum, violating one of 
the key assumptions made in calculating MBondi- This is 
important because it implies that infalling material will 
settle onto a circular orbit whose radius i? c irc is set by the 
angular mome ntum of the mate rial with respect to the 
black hole (cf. iHobbs et~aT1 l2010h . In particular, it means 
that only the very lowest angular material will be available 
to feed the black hole because the timescale required for 
viscous transport of material through the disc is of order 
a Hubble time on sc ales of order R ~ 1 — lOpc (see, for 
example, lKindl2010h . This is a very restrictive condition 
because it is not straightforward for infalling gas to lose its 
angular momentum other than by colliding with other gas, 
which leads to angular momentum cancellation. Therefore, 
angular momentum provides an efficient natural barrier to 
accretion by the black hole, and so must be accounted for 
when estimating Mbh- 

These arguments make clear that the Bondi-Hoyle 
method cannot provide a reliable estimate of Mbh in galaxy 
formation simulations. If feedback from black holes plays 
as i mportant a role in galaxy formation a s we expect it to 
fe.g. lBower eTa l. 2006; Croton et al.ll2006l \ then it is crucial 
that we devise an alternative method for estimating Mbh in 
galaxy formation simulations that overcomes the problems 
that beset the Bondi-Hoyle method. 

In this short paper, we introduce our new "accretion 
disc particle" method (hereafter the ADP method) for es- 
timating Mbh in galaxy formation simulations, which ac- 
counts naturally for the angular momentum of infalling ma- 
terial. We use a collisionless accretion disc particle (ADP) 
to model the black hole and its accretion disc. The black 
hole accretes if and only if gas comes within the accretion 
radius R acc of the ADP, at which point it is captured and 
added to the accretion disc that feeds the black hole on a 
viscous timescale t v i SC - In this way the black hole will ac- 
crete only the lowest angular momentum material from its 
surroundings. 

The layout of this paper is as follows. We describe the 
main features of the ADP method in £[2] showing how the ac- 
cretion rate M acc onto the ADP is linked to the black hole ac- 
cretion rate Mbh - In fc|2. 2 1 we discuss briefly our mom entum- 
driven feedback model (cf . iNavakshin &: PowerfeplOh as well 
as our implementati on of the quasar pre-heating model of 
ISazonov et al.l ([2005T ) . The accretion rate Mbh estimated us- 
ing the ADP method is very different from one estimated us- 
ing the Bondi-Hoyle method. We show this clearly in fus- 
ing simple idealised numerical simulations, designed to illus- 




SPH Particle B 



Figure 1. The Accretion Disc Particle (ADP) Method. 

The ADP is a collisionless sink particle that consists of a black 
hole and its accretion disc. SPH particle A has a small angular 
momentum and so its orbit brings it within the ADP's accretion 
radius R aC c, at which point it is added to the accretion disc. SPH 
particle B's angular momentum is too large for it to be captured. 
The black hole feeds from the accretion disc on a viscous timescale 
t v i sc . Both i^acc and £ v i sc are free parameters in the accretion disc 
particle method. See text for further details. 

trate the key differences between the ADP and Bondi-Hoyle 
methods for estimating Mbh- These simulations follow the 
collapse of an initially rotating shell of gas onto a black hole 
embedded in an isothermal galactic potential. Finally we 
summarise our results in gland we discuss the implications 
for galaxy formation simulations and the Mbh — cr relation 
in 23 

2 MODELLING ACCRETION AND 
FEEDBACK 

2.1 The Accretion Model 

2.1.1 The Accretion Disc Particle (ADP) Method 

The main features of the ADP are illustrated in Fig [1] 
The ADP is col lisionless and consists of a sink particle 
(jBate et al.lll995h with an accretion radius R acc • Race is a 
free parameter of the simulation but in general it is desirable 
to set it to the smallest resolvable scale in the simulation, 
which will be of order the gravitational softening length of 
the gas particles. The total mass of the sink particle is equal 
to the sum of the masses of the black hole, Mbh, and its 
accretion disc, Mdi sc - The accretion disc is assumed to be 
tightly bound to the black hole and is thus a property of the 
sink particle rather than a separate entity. 

Accretion onto the black hole in the ADP method is a 
two- stage process. First, any gas that crosses the accretion 
radius i^acc is removed from the computational domain and 
added to t he accretion disc. In the classical sink-particle 
method of iBate et all (|l995h , the accreted gas would be 
added to the black hole immediately, but in an astrophysi- 
cally realistic situation, the finite non-zero angular momen- 
tum of the accreted gas leads to the formation of an accretion 
disc before the gas can accrete onto the black hole. Here we 
assume that gas is added to the accretion disc after a time 
that is of order the dynamical timescale tdyn at R acc . 

Second, gas is transported through the accretion disc 
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and is added to the black hole. In principle, we could de- 
scribe the evolution of the accretion disc by the standard 
viscous disc evolution equa tions (see, for example, Chap- 
ter 5 of iFrank et al.1 l2002h . However, neither theory nor 
observation tell us what the magnitude of the disc vis- 
cosity should be and so we cannot be sure of the effi- 
ciency of angular momentum transport within the disc. 
Moreover, if the accretion disc is sufficiently massive to 
become self- gravitating, then a significant fraction of its 
gas mass can be converted to stars (see, for example , 
Toomrd[l964l: |Paczvhskilll978l ; IShlosman fc Begelman 1989; 



Goodman! I2OO3I : iNavakshin et alJl2007h . Star formation de- 
pletes the gas reservoir, reducing Mbh, but st ellar feedback 
can act to either reduce or en hance Mbh (e.g. ICuadra et al.l 
120081 : ISchartmann et alJl2009h . Solving the standard viscous 
disc evolution equations is no longer straightforward under 
these circumstances and detail ed accretion disc sim ulations 
are required (see, for example ICuadra et alj |2008). 

We want a simple sub-grid model, however, and so we 
simplify the problem. To this end, we assume that (a) an 
accretion disc forms and (b) angular momentum transport 
through the disc introduces a delay between the time a gas 
particle crosses R acc and the time that it is accreted by the 
black hole. This time delay will be of the order of the disc 
viscous time, t v isc, which can be of order the Hubble ti me for 
accre tion discs around super-massive black holes (cf. iKind 
To capture this in the simplest way we describe the 
evolution of the accretion disc by 



M disc = M acc - M BH , 



(1) 



where M acc is the rate at which gas is captured through R acc 
and Mbh is the accretion rate onto the black hole. t v i sc can 
be es timated using physical arguments (see, e.g. §2 of iKind 
l201Ch but, as argued above, it is uncertain, and so we treat 
tvisc as a free parameter instead. In general we require that 



tdyn(^acc) and typically set t v 



10 - lOOMyrs in 



tv 

galaxy formation simulations (Power et al., in preparation). 

Note that the rate at which gas is captured from the 
ambient medium is not limited in any way; it is simply gov- 
erned by the evolution of the large scale accretion flow into 
the centre of the galaxy. In contrast, the rate at which gas 
is accreted onto the black hole is limited by the Eddington 
accretion rate MecM, 

• _ 4tt G Mbh m p 

M E dd = , (2 J 

?7 (XT C 

where m p is the proton mass, <tt is the Thomson cross- 
section, c is the speed of light and r\ is the accretion effi- 
ciency, for which we assume the standard value of r\ = 0.1. 
This means that Mbh in equation [1] satisfies 

"Mdisc 



Mbh 



tvis 



,M E 



(3) 



This simple system of equations can be expanded in the 
future to encompass more detailed disc modelling, including 
gas disc self-gravity and the resulting star formation and 
feedback from stars formed there. Gi ven the empirical evi- 
dence from our Galactic C entre (e.g., INavakshin &: Cuadral 
120051 ; [Paumard et al.ll20 06) and the theoretical expectation 
tha t nuclear stellar cluste r feedback should be important 
(cf . INavakshin et alJl2009bh , this is likely to be an important 
step in future studies. 



It is worth making some additional comments about our 
ADP method and how it relates to the classical sink particle 
formulation o f iBate et al.1 (jl995T ). In the classical sink parti- 
cle formulation, a number of conditions had to be satisfied 
before gas could be accreted by the sink particle (e.g. pres- 
sure forces at the accretion radius, comparison of thermal 
and gravitational binding energy with respect to the sink 
particle, etc.). In our approach, there is a single condition 
for accretion, namely that gas comes within R acc - Physically 
this is quite reasonable. The virial temperature in the vicin- 
ity of a super-massive black hole is very high, typically in 
the range of 10 6 — 10 8 K. However, gas densities are also very 
high near R acc and so cooling times are expected to be very 
short (e.g. King 20051 ). This implies that gas is likely to be 
much cooler than the virial temperature, which means that 
both the pressure forces and thermal energy of the gas is 
negligible. Furthermore, viscous times are always very long 
compared with dynamical times, so we expect the accre- 
tion disc to be a long-lived (essentially permanent) feature 
within R a cc- This means that gas that comes within R acc is 
very likely to undergo a large Mach number collision with 
the disc, causing it to shock and then cool rapidly. There- 
fore, even if gas is initially on an unbound (hyperbolic) orbit 
around the sink particle, it will most likely lose most of its 
bulk and thermal energy and settle into the accretion disc. 



2.1.2 The Bondi-Hoyle Method 

The Bondi-Hoyle method for estimating Mbh is the 
standard approa ch in galaxy formation s imulation s (see , 



for example, iDi Matteo et al.1 120051 : ISpringell [2005; 



iBooth fe S chave 20091). Here the black hole accretion rate 



is calculated directly from 



Mb 



47raG 2 Mj H p 

( c 2 +v 2)3/2 ' 



(4) 



where p is the SPH density at the position of the black hole, 
c s is the sound speed of the gas, v is the velocity of the black 
hole relative to the gas and a is a fudge factor that we set to 
unity for the purposes of this work, but which can be of orde r 
- 100 - 300 (see the discussion in IBooth fe Schavd l20Q9h . 
In practice we compute estimates for p, c s and v using the 
SPH smoothing kernel with A/sph=40 neighbours. Note that 
there is no explicit dependence on the angular momentum 
of the gas in equation 2] - the accretion rate is dictated by 
the gas density p and sound speed c s . 



2.2 The Feedback Model 

In the simulations presented in the next section, we use Mbh 
estimated using either equation [3] or U to determine the ac- 
cretion luminosity of the black hole, 



r]M B uc 



(5) 



this is Eddington limited, as explained in the previous sec- 
tion. We assume that this radiated luminosity drives a wind 
that carries a momentum flux L acc /c, which is usually true 
for AGN (cf. iKing fe Pounds] 120031 : lKindl2010h . Wind par- 
ticles are emitted isotropically by the black hole at a rate 



-N wind 



CPwind 



(6) 



© 2010 RAS, MNRAS 000, EH 



4 C. Power, S. Nayakshin & A. R. King 



and they carry a momentum p w ind = O.lra ga s0", where a is 
the velocity dispersion of the host halo. This satisfies the 
requirement p w ind <C p gas , where p gas is the typical gas par- 
ticle momentum (~ 

eras 

cr here), and ensures that Poisson 
noise from our Monte Carlo scheme d oes not compromise 
our results fsee iNavakshin et al.ll2009al ). 

In addition to this momentu m-driven wind, we inc lude 
the quasar pre- heating model of ISazonov et all (|2005l ). In 
this model , the average quasar spectral energy distribution 
derived by Sa zonov et al is used to estimate an equi- 

librium temperature T eq for the gas based on the ionisation 
parameter £(r) = L acc /n(r)r 2 , where n(r) is the number 
density at radius r. Physically T eq corresponds to the tem- 
perature at which heating through Compton scattering and 
photoionisation balances Compton cooling and cooling as a 
result of continuum and line emission, on the assumption 
that gas is in ionisation equilibrium. In practice, we calcu- 
late heating and cooling; rates using fo rmulae A32 to A39 in 
Appendix 3.3 of I Sazonov et al.l (|2005h . and we find that the 
resulting equilibrium temperature profile of the gas is well 
approximated by their equation 3, 

T eq (0 - 200£K. (7) 

This holds over the temperature range 2 x 10 4 K to 10 7 K; 
for f < 100 and f > 5 x 10 4 , T eq ~ 10 4 K and 2 x 10 7 K 
respectively. 



3 RESULTS 

We have run simple idealised numerical simulations that are 
designed to show that our ADP method constitutes a physi- 
cally self-consistent sub-grid model for estimating the black 
hole accretion rate Mbh in galaxy formation simulations, 
and that the Bondi-Hoyle method does not. 

Our initial condition is a spherical shell of gas of 
a uniform density po, distributed between the inner and 
outer radii, Ri n and R ou t, respectively. The shell is em- 
bedded in the static analytic gravitational potential of a 
singular isothermal sphere with a 1-D velocity dispersion 
a and modified slightly to have a constant density core 
within R ^ R C ore- For all the runs in this paper, we adopt 
po ~ 10 10 M©kpc~ 3 , R in = 0.067 kpc, R ou t = 0.1 kpc, 
^Rcore = 0.01 kpc and a — 147kms _1 . The shell has a mass 
of M sh eii = 3 x 10 7 M o and is realised with ~ 280, 000 gas 
particles, drawn from a uniform density glass, which means 
that the particle maSS IS TTlp-as ~ 1.1 x 10 2 M . We give the 
shell an initial temperature of 10 4 K and an initial bulk ro- 
tation around the z-axis such that its rotational velocity in 
the x-y plane is v<f> — v ro t — /rot V%cr with / ro t = 0.3; it falls 
from rest in the radial direction. Finally we embed a colli- 
sionless particle - corresponding to the black hole - at rest 
at the centre of the potential; the initial black hole mass is 
Mbh = 10 6 Mq and, in the cases where we use the accretion 
disc particle model, an initially zero disc mass. 

All of the simulations are run us ing GADGET3, an up- 
dated version of the code presented in lSpringel (|2005h . Each 
simulation is run for ~ 4.7 Myrs, which corresponds to ~ 14 
dynamical times at the initial outer radius of the shell. 




-4-2 2 

log E [g cm" 2 ] 




-4-2 2 

log E [g cm" 2 ] 



Figure 2. Without Feedback : The projected gas density in 
the x-y and x-z planes at £=4.7 Myrs (upper and lower panels 
respectively). As before, arrows indicate the velocity vectors of the 
gas. Note the dense knots corresponding to high density regions 
of sink particle formation. The ADP has an accretion radius of 
i?acc = 0.003 kpc. 



3.1 Without Feedback 

We begin by considering the simplest possible case - the 
collapse of the shell in the absence of any feedback from the 
black hole. We model the black hole using an ADP, but in 
this particular simulation we decouple the accretion lumi- 
nosity Lace of the black hole from Mbh by setting r\ = 
in equation [5] this suppresses both the momentum-driven 
wind and quasar pre-heating. We choose an accretion radius 
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of R a cc = 0.003 kpcQ- For simplicity, we assume an isother- 
mal equation of state with a temperature of T = 10 4 K. 

By conservation of angular momentum, th e shell should 
settle into a thin rotationally supported disc (cf. lHobbs et al.l 
l20ld ). This disc is distinct from, and on a much larger scale 
than, the accretion disc discussed in ^2.11 which is tightly 
bound to the black hole on a scale much smaller than we can 
resolve in our simulation. We show the gas density projected 
onto the x-y and x-z planes at £=4.7 Myrs (upper and lower 
panels respectively) in Fig [2] The gas is distributed in a thin 
rotating disc; its inner and outer boundaries are at ~ 0.006 
and ~ 0.01 kpc respectively and it rotates in a clock-wise 
sense around the z-axis (indicated by the projected velocity 
vectors). The inner boundary is larger than the accretion 
radius R acc = 0.003 kpc by a factor of ~2 and so only a 
small fraction of the mass of the disc comes within R 

acc over 

the duration of the simulation (~20 particles or ~ 0.007% 
after ~ 4 Myrs). 

The absence of any accretion until late times might seem 
counterintuitive, given that our initial condition is a rotating 
shell of gas. What happens to material with small angular 
momentum that lies along the axis of rotation? Why d oes it 
not accrete rapidly onto the ADP? lHobbs et al.l (|2010h have 
shown that this small angular momentum gas shocks and 
mixes with larger angular momentum gas, which increases 
its net angular momentum and provides a barrier to accre- 
tion. 

Note that there are knots of high density material in 
the disc. For expediency, we tag gas particles that exceed a 
threshold value in local density. Physically these high density 
regions are likely to host star formation, but for the purpose 
of this study we simply decouple these particles hydrody- 
namically from other gas particles, ignoring them in the hy- 
drodynamical force calculation, which helps to increase the 
speed of the simulation. In this particular run, ~ 90% of the 
gas particles are converted to decoupled particles by t ~ 4.7 
Myrs. 

3.2 With Feedback 

Let us now consider the evolution of the collapsing shell 
when the black hole accretion luminosity L acc is coupled to 
Mbh, as in equation [5] We estimate Mbh using the ADP 
method (equations [T] and [3] with R acc = 0.003 kpc and 
t v i sc ~ 10 4 yrs) and the Bondi-Hoyle method (equation [H 
with a=l). Note that, for the purpose of this study, the pre- 
cise value of tvisc is unimportant; the point is that in the 
cases that we consider, the ADP method predicts a negligi- 
ble accretion rate, as we would expect from physical argu- 
ments. In both cases we assume that the black hole feedback 
takes the form of a momentum-driven wind and quasar pre- 
heating. The gas is initially isothermal with a temperature 
of T = 10 4 K, but as it evolves it can heat and cool in 
response to, for example, the quasar radiation field. 

In Fig. [3] and Fig. [4] we show the gas density projected 



1 This is larger than the gravitational softening of gas particles 
by a factor of ~ 20, but slightly smaller than the inner edge of 
large-scale gas disc we expect to form. We choose a large value to 
highlight that our initial conditions minimal accretion onto the 
black hole, by construction. 



onto the x-y and x-z planes (upper and lower panels respec- 
tively) in the ADP and Bondi-Hoyle runs (left and right 
panels) at t—1 Myr and £=4.7 Myrs. As in Fig. [2j arrows 
indicate the magnitude and direction of the projected ve- 
locity vectors of the gas. The differences between the runs 
are striking. The shell should settle into a thin rotationally 
supported disc whose properties are very similar to those of 
the disc shown in Fig. [2] and indeed this is the case in the 
ADP run. This is unsurprising - because we link feedback 
explicitly to accretion rate onto the black hole, we do not 
expect any significant feedback in the ADP run because the 
angular momentum of the gas is too large to bring it within 
Race until late times. Even at this point, the mass of gas 
accreted M acc <C 1% over the lifetime of the simulation and 
Mbh <C MecM, which means that the feedback is weak and 
has little effect on the gas distribution. We note also that, 
as in the run without feedback, ~ 92% of the gas particles 
have decoupled into high density knots by t c± 4.7 Myrs. 

In contrast, the accretion rate is consistently Eddington 
limited over the duration of the simulation in the Bondi- 
Hoyle run. The black hole is not massive enough for its feed- 
back to prevent the collapse of the shell into a disc (cf. right 
hand panels of Fig. [3]), but once the disc has formed, the 
feedback acts efficiently on the low column density gas. At 
early times it is the lower column density gas surrounding 
the disc and along the axis of rotation that is most efficiently 
driven outwards, principally by the momentum-driven wind. 
In particular, it is the impact of the feedback along the axis 
of rotation that imprints the strongly bipolar character on 
the outflow, evident in the x — z projection in Fig. [3] Over 
time, as high density knots form in the disc, lower density 
material within the disc is blown away. For example, after 
t ~ 1/2.8/4.7 Myrs, -2%/25%/40% of the gas particles that 
have not been decoupled have been expelled from the disc, 
compared to <C 1% of gas particles in the accretion disc 
particle run. At late times, even the high density knots are 
ablated. As a result, the cumulative effect of the feedback 
over ^4.7 Myrs has a profound impact on the gas distribu- 
tion, a point that is nicely illustrated in Fig.[5l which makes 
clear that the gas is distributed over kpc-scales (and indeed 
to ~ 10 kpc) by the end of the run. 

It is worth noting that the fraction of gas that had a 
sufficiently high density to be decoupled over the lifetime of 
the Bondi-Hoyle run is comparable to the fraction in the ac- 
cretion disc particle run (88% compared to 90%). However, 
these high density knots are ablated by the feedback in the 
Bondi-Hoyle run whereas they survive in the accretion disc 
particle run. 



4 SUMMARY 

Black holes grow by accreting gas and stars from their sur- 
roundings. However, only the lowest angular momentum ma- 
terial can come sufficiently close to the black hole to be ac- 
creted, and so any estimate of a black hole's accretion rate 
Mbh must account for this. However, the standard sub-grid 
model for black hole accretion in galaxy formation simula- 
tio ns neglects the angular momentum of accre ting material 
fcf. lDi Matteo et a l. 2005; Springel e t al. 2005). The Bondi- 
Hoyle method (Bondi Hovld 19441 ; lBondilll952h assumes 
that Mbh oc pj(? s where p is the gas density at the position 
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Figure 3. Early Times: The gas density projected onto the x-y and x-z planes (upper and lower panels) in the accretion disc particle 
and Bondi-Hoyle runs (left and right panels) at t ~ 1 Myr. 



of the black hole and c s is the sound speed in the gas. This 
implies that black holes are always accreting; Mbh may be 
small but it can never be zero, regardless of the angular 
momentum of the gas surrounding the black hole. 

In this short paper, we presented a new sub-grid model 
for estimating Mbh in galaxy formation simulations that 
accounts for the angular momentum of accreting material. 
This accretion disc particle (ADP) model uses a collisionless 
sink particle to model the composite black hole and accre- 



tion disc system. The black hole accretes if and only if gas 
comes within the accretion radius i^acc of the ADO, at which 
point its mass is added to the accretion disc that feeds the 
black hole on a viscous timescale t v i sc - In this way the black 
hole will accrete only the lowest angular momentum mate- 
rial available to it in and around its host galaxy. 

We demonstrated that the ADP method constitutes a 
physically self-consistent model using simple idealised nu- 
merical simulations that follow the collapse of a rotating 
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Figure 4. Late Times: The gas density projected onto the x-y and x-z planes (upper and lower panels) in the accretion disc particle 
and Bondi-Hoyle runs (left and right panels) at t ~ 4.7 Myrs. 



shell of gas onto a black hole embedded at the centre of 
an isothermal galactic potential. By construction, the gas 
settles into a thin rotationally supported disc between Ri n 
and i?out, where R m in > Race, and so we do not expect any 
accretion onto the black hole. Because we link feedback to 
accretion, we do not expect there to be any significant dif- 
ferences between simulations with or without feedback when 
using the ADP estimate of Mbh- On the other hand, we ex- 



pect the evolution of the system to differ if Mbh is estimated 
using the Bondi-Hoyle method. 

These expectations were borne out by the results of our 
simulations. The Bondi-Hoyle method predicted that Mbh 
should be Eddington limited over the lifetime of the simula- 
tion. Because the feedback in this case was relatively weak, 
its effect could not prevent the collapse of the shell into a 
disc, but the cumulative effect of the feedback was to drive 
gas away and to expel it to ~ lOkpc scales after ~ 5 Myrs. 
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Figure 5. Large-scale distribution of gas in the Bondi-Hoyle run 
at t ~ 4.7 Myrs. 



In contrast the ADP method predicted negligible accretion 
rates at all times; the shell collapsed, settled into a thin 
rotationally supported disc and ~ 90% of the mass decou- 
ples into long-lived high density knots, which correspond to 
regions of star formation. 



5 CONCLUSIONS 



We have argued that our new accretion disc particle (ADP) 
method provides a far more physically motivated and 
self-consistent approach to modelling black hole accretion 
than the Bondi-Hoyle method, which is th e standard ap- 
proach in galaxy formation s imulations (cf. ISpringel et al.l 
120051 : iDi Matteo et all l2005h . The Bondi-Hoyle method 
was formulated with a specific astrophysical problem in 
mind, quite unlike the situations that arise when modelling 
galaxy formation. It is not applicable to problems in 
which the accretion now has non-zero angular momentum 
(as demonstrated in this paper) and/or in which it is 
embedded in the potential of a more massive host (as 
we show in Hobbs et al., in preparation). Therefore it is 
unsurprising that the Bondi-Hoyle method struggles to 
capture the behaviour of gas accretion in these kinds of 
common situations. Our ADP method is similar in spirit 
to the "accretion radius" or "sink particle" approaches to 
modelling accretion that are used extensively in simulating 
star formation (e.g., iBate et al.lll995l : feate fc Bonnellll2005h 
and modelling gas accretion onto the super-massive blac k 
hole at the centre of the Milky Way (|Cuadra et alJl2Q06h , 
and we believe that it is natural to extend this approach 
into modelling galaxy formation. An important next step 
in our work, which builds on this accretion disc particle 
method and our recent momentum-dri ven wind model for 
feedback (cf. iNavakshin fe Powerll2010h . is to combine the 
models in simulations of merging galaxies and ultimately 
cosmological galaxy formation simulations. 



It is interesting to consider one important astrophysical 
consequence of our accretion disc particle model and to con- 
trast it with what one would expect using the Bondi-Hoyle 
model. Recently it has been suggested that there is observa- 
tionally evidence for separate Mbh — cr relations for elliptical 
galaxies and classical bulges on the one hand and pseudo- 
bulges on the other, such that the bla ck holes in pseudo- 
bulges ar e underweight (see, for example jGreene et al 1 120081 : 
iHul [2009). The properties of pseudo-bulges appear to de- 
viate systematically from those of classical bulges, and in 
particular they are characterised by a high degree of rota- 
tional support. As we have shown, the angular momentum 
of infalling material provides a natural barrier to black hole 
growth, and so we would expect that rotationally supported 
systems to be more likely to be systems in which the central 
super-massive black hole is malnourished and underweight. 
Of course, the precise details of a galaxy's assembly history 
are important but our model would predict a systematic off- 
set between super-massive black hole masses in galaxies that 
have on average accreted higher angular momentum mate- 
rial than those that have on average accreted lower angular 
momentum material. In contrast, the Bondi-Hoyle model 
would predict that the black hole should continue to grow 
to a critical black hole mass imposed by the depth of the 
gravitational potential in which it sits, regardless of the an- 
gular momentum of infalling material. We shall investigate 
this question further in future work. 
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